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Abstract 

Exact eigenvalue correlation functions are computed for large N her- 
mitian one-matrix models with eigenvalues distributed in two symmetric 
cuts. An asymptotic form for orthogonal polynomials for arbitrary polyno- 
mial potentials that support a Z2 symmetric distribution is obtained. This 
results in an exact explicit expression for the kernel at large N which de- 
termines all eigenvalue correlators. The oscillating and smooth parts of the 
two point correlator are extracted and the universality of local fine grained 
and smoothed global correlators is established. 



1 Introduction 



Eigenvalue correlators in large random matrix models have been of interest 
recently (see [Q, [|, f|, [|, ^, |7], ^, ^ and references therein). The recent 
interest is in part inspired by the burst of activity in the field of quantum 
chaos and mesoscopic systems. 

Correlators in matrix models have been calculated in various limits by 
different methods. At short separations of the order of the eigenvalue spacing 
the correlators show oscillatory behaviour. Correlators that reproduce this 
behaviour are referred to as "fine grained" . If these oscillations are averaged 
over, one gets the so called "smoothed" correlators. For some applications, 
one needs the correlation functions over the entire eigenvalue distribution, 
near the centre as well as the edges (global correlators) while for others only 
in a small region near the centre or the edge of the eigenvalue distribution 
(local correlators). The complete information is contained in the global fine 
grained correlator, from which all others can be obtained as special limiting 
cases. 

The global fine grained correlator has been computed for the hermitian 
one matrix model []!]] and the hermitian two matrix model [10] using the 
method of orthogonal polynomials. These studies, however, consider only 
those cases where the eigenvalue distribution lies in a single cut. The present 
paper is concerned with correlation functions for hermitian matrix models 
in which eigenvalues lie in two cuts or bands. 

Random matrix models with eigenvalues that have support in more 
than one cut have been studied in the context of 2-d quantum gravity and 
symmetry breaking solutions (pT[ and references therein), two-dimensional 



QCDjT^, [L3| and statistical physics They may also be applicable in 

studies of certain condensed matter systems with eigenvalue distributions 
with gaps. In the present paper the global fine grained correlators for matrix 
models with eigenvalue support in two symmetric cuts are obtained. These 
expressions therefore contain detailed information about how the correlators 
vary over distances of the eigenvalue spacing, from which the smoothed cor- 
relators can be extracted by a suitable averaging procedure. They are also 
valid over the entire eigenvalue distribution, near the centre as well as the 
edge of the cuts. 

In order to compute the global fine grained correlators, an asymptotic 
form for orthogonal polynomials is derived, which is valid for arbitrary po- 
tentials that support a Z2 symmetric eigenvalue distribution. This form 
allows a determination of an explicit expression for the eigenvalue kernel 
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through the Christoffel-Darboux formula suitably modified for the two cut 
situation. The kernel then determines all the global fine grained n-point 
correlation functions of the eigenvalue density. 

The two point function is explicitly discussed in various limits. The 
kernel as well as the global fine grained correlator in general depend upon 
the coupling constants in the matrix model potential. However, it is shown 
that in two limits the two point function is universal, i.e., independent of the 
form of the potential. These two limits are (i) fine grained local correlators 
away from the edge of the cuts, where the Dyson form is reproduced, and 
(ii) global smoothed correlators. 

This generalizes the work of to the two cut matrix models. When 
the coupling constants of the matrix model change such that the eigenvalue 
support moves from two cuts to one cut, the single-cut results are obtained. 
Thus the present work can be used to study this phase transition in the cou- 
pling constant space of matrix models at the level of fine grained correlators, 



generalizing earlier studies at the level of smoothed correlators [jD], 16]. 

The paper is organized as follows. In the second section after estab- 
lishing notations and conventions, the asymptotic form for the orthogonal 
polynomials for two symmetric cuts is presented as an ansatz. In the third 
and fourth sections the recurrence relation and the orthogonality condition 
are used to determine the unknown functions in the ansatz. The fifth sec- 
tion derives the Christoffel-Darboux formula for the kernel K(fi, v) in a 
form suitable for two cuts and an explicit expression for kernel is presented. 
The sixth section discusses the two point correlation function in its full fine 
grained global form and establishes its universality in various limits (fine 
grained local, smoothed global). The seventh section contains concluding 
remarks. 



2 Notation, Conventions and Asymptotic Ansatz 

The partition function is defined by 

Z = J dM exp(-iV TrV(M)) , (2.1) 

where M is an N x N hermitian matrix and V(M) = Ysr£=i ^M n is a 
real polynomial in M (the same notation as ref. JTTj ] is used). P n (X) is a 
set of polynomials that are orthogonal with respect to the measure defined 
by the potential V, J^° 00 dXP n (X)P m (X)exp(—N TrV(A)) = h n 5 nm , and are 
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normalized such that P n (X) = A n +ciA n_1 + ...., Po(X) = 1. These properties 
determine P n (X) uniquely for every V. The orthogonal polynomials satisfy 
the recursion relation AP„(A) = P ra+ i(A) + S n P n (X) + R n P n -i(X) where the 
recursion coefficients R n , S n are V dependent. R n = h n /h n -\ and for Z2- 
symmetric potentials (V(— A) = V(A)) S n = 0. It is convenient to define an 
orthonormal set of polynomials ip n (\) = Pn 2^ exp{— ^ Try (A)). Then 

J dX^ n (X)^ m (X) = 5 nm , (2.2) 

and 



A^(A) = /R^Vn+i(A) + V^V'n-i(A). (2.3) 
The kernel K(p, v) is defined as 

-. N-l 

K(v,v) = -J2^(P)M")- (2-4) 

i=0 

The Christoffel-Darboux formula expresses this in terms of just ip n with n 
close to N: 

K(H, v) = • (2.5) 

The kernel determines all eigenvalue correlators. Defining 
p{p) = jfTrS(p ~ M) = jfJ2iLi^{p ~~ ^i)) where Aj are the eigenvalues 
of M, it follows that 

p{p) = (p(p)) = k (p,p) 
Pc {li,v) = (p(p)p(v))-(p(p)){p(v)) = -[K(p,v)} 2 

p c {\,H,v) = {p(X)p(p)p(v)) c = 2K(X,p)K(p,v)K(v,X) (2.6) 

where the subscript c stands for connected correlator and the expectation 
value of any function f(M) is given by (/) = \ J dMe~ NTrV ^ f(M). 

When the potential is Z2 symmetric and has two wells the eigenvalue 
density p(X) can have support in two cuts located symmetrically about the 
origin, A € (a, b)U(—b, —a), a > 0, b > a. One is then interested in the kernel 
and correlation functions where p, u, etc. lie in the cuts. The Christoffel- 
Darboux formula implies that we need the orthogonal polynomials for n 
close to N. 
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For A lying in the two cuts and for N — n ~ 0(1), N large, we make the 
ansatz that the orthogonal polynomials can be approximated by 

^„(A) = -^j[cos(NC -(N-n) ( / ) + x + (-lfr?) + O(l)], (2.7) 

where /, £, (j>, X an d r\ are functions of A, and that tp n is damped out outside 
the cuts. We then show that 

/(A) = |^^^ sin 20(A) (2.8) 

from the orthogonality condition satisfied by the orthogonal polynomials 
and is derived in section |||. £(A) is fixed to be 

C'(A) = -ttp(X) (2.9) 

from the relation K(X, A) = p(X) and is derived in section [|. 0(A) and 
77(A) are determined from the recurrence relation satisfied by the orthogonal 
polynomials ( section [| ) 

A 2 _ (a 2 +b 2 ) 

cos20(A) = — , (2.10) 
2 



, cos 0(A) 

cos2r?(A) = b p~i, 

A 

sin (6(A) 

sin 217(A) = a ^-A (2.11) 

A 

This leaves the function x(^) undetermined and we conjecture that it has 
the form — f for all potentials as in the single cut case ref. ([jl], |l7fl). 

Example for #2 < 0,<?4 > 0, all other gi = 0, and in the two cut phase 
g| > 4^4, the end points are given by a 2 = — 2^/J — ^ and 6 2 = 2y^- — 
The density of eigenvalues is p(A) = |^A V / (A 2 — a 2 )(b 2 — A 2 ) for A G 

(—6, —a) U (a, 6) and zero elsewhere. Hence ((X) = 0(A) — sin4 ^ A ) t This 
determines the asymptotic form of ip n (\) completely. 

Note that /, 0, 7] and x are universal functions, i.e, their functional form 
is independent of the potential V, the only dependence on V enters through 
the endpoints of the cuts a and b. £ is non-universal since the eigenvalue 
density depends in general on the detailed form of V. 
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3 Recurrence Relation 



We will now use Eq. (|2.3| ) to determine 4> and r\. Multiplying Eq. (| 
A and using Eq. (|2.3D once again we get 



A 2 V'n(A) = v 7 Rn+1 V Rn+2lpn+2 (A) + R n+ l1p n (X) + R n ljj n (X) 



+ 



Rn V -Rn-lYVi-2 ( A) • 



by 



(3.1) 



When V has two symmetric wells, it is known that Rn = A n for n even and 
R n = B n for n odd, where in the large N limit A n and B n are approximated 
by two continuous functions A(x), B(x) with x = -8= [18, |l(|. Thus, for n 
even 



A 2 ^n(A) = y/B n+ i ^ A n+2 lfj n +2(X) + B n+1 1p n (\) 
+ A n i> n (\) + V / Xv / ^n-l^n-2(A). 



(3.2) 



On using the asymptotic ansatz for even n i.e., substituting Eq. (|2.7| ) in Eq. 
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JV-ji 



(|3.2|), and replacing A n = A(x = -8>) = A(x — ^ — ^ 
A(l) = A (and similarly ,B n — > 5(1) = 5) we get 



)=A(l)+0(4 



A 2 - (A + 5) 



cos 20(A). 



(3.3) 



2VAB 

It is known that for Z 2 symmetric two cut models the end points a and b 
are related to A and B by [16] 

A + B 

2v / AB 



a 2 + 6 2 



(3.4) 



This yields Eq. ( pTlOj) . 

Next consider the recurrence relations (n even) 



A^ n+ i(A) 

A^n-l(A) 



V ARn+2'0n+2(A) + v / -Rn+l'0n(A) 
y/R^1p n (X) + v / -Rn-lV'n-2(A). 



(3.5) 

On substituting Eq. ( |2.7| ) for ?/;„, ip n ±i, ip n ±2 in Eq. ( |3.5D it is easy to see 
that 



sin 277(A) 



(-Va + Vb) 

sin 0(A) 



sin 4>(X) 
A 



(3.6) 



5 



and 

cos 277(A) = (VA + Vb)^^ 

A 

^cos (p(X) 
A 

This determines Eqs. ( [2.10D and (2.11). 



(3.7) 



4 Orthogonality 



Let us check for orthonormality, Eq. (2^2), for all n,m. The ansatz assumes 



that for large ./V t[i n (\) is damped out sufficiently fast beyond the end of 
cuts so that integration can be restricted to the cuts 



( / + j h )d\^ n {X)^ m {\) = 5 nm . (4.1) 



The integrand is 



1 

WnVm = J COS a COS 

= — [cos(a + 13)+ cos(a - 0)] (4.2) 
2/ 

where (a + 0) = (2N( - (2N - n - m)<f> + 2 X + (-1)» and (a - 0) = 
{{n — m)(p+{{— l) n — (— l) m )rj). The cos(a + 0) term is zero upon integration 
because cos(2iV£ — (2N — n — m)4> + 2\ + (— l) n r/) oscillates rapidly and 
averages to zero. Thus 

b + J a ) T cos{a-0) (4.3) 

Note that (-l) n - (-l) m = if n - m is even and (-l) n - (-l) m = 2(-l) n 
if n — m is odd. For n — m even 

f—a rb ^ 

(V'n^m) = (/ + / )— cos(n-m)0. (4.4) 

J-b Ja 2/ 

In particular for n — m = 

/—a rb ,7\ 

-b + Ja ^ (45) 
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We will show that if 

37 = C /( A )> ( 4 - 6 ) 



then orthogonality follows. On taking the derivative of Eq. ( 2.1C| ) one gets 

dX -1 {b 2 - a 2 ) 



d(j> A 

Thus 



■sin 2^ (4.7) 



f W = Z* (*-«'>*■*(*), (4.8) 
c 2 A 

Note from Eq. (|4.5| ) that for (ip n ,if) n ) = 1, /(A) must be an even function of 
A. This implies that sin 20 must be an odd function of A. Since cos 20 = —1 
at A = ±a, and 1 at A = ±6, this is achieved by taking 20 = 2pir, (2p — 
l)"7r, (2p — 3)7r, (2p — 4)tt at A = —b, —a, a, b respectively; with p any integer. 
With this choice sin 20 is odd because 20 is in the first two quadrants for A 
in (a, b) and the third and fourth quadrants for A in (—6, —a). In particular 
this choice also implies that cos is an even function of A and sin is an 
odd function. Substituting Eq. ( |4.8| ) into Eq. Q4.5D implies that c = — . 
Using c = ^ in Eq. ( jig ) gives Eq. (gj). 



Now consider n — m ^ even. Using Eq. (|4.6|) in Eq. (4.4), we find 



r<t>{-a) r<t>{b) 
(ipniipm) = c( / +/ )#cos(n - m)0 

J</>(-6) J<Ma) 

= 0. (4.9) 

Though the integrand is an even function of A the full integral vanishes as 
each integral separately vanishes. 
Now consider n — m odd 

+ / )— cos(a-/?) (4.10) 

-b Ja 2} 

where a — [5 = (n — m)0 + 2(— l) n r/ for n — m odd. Expanding cos(a — (3) 
we get 

r-a rb d\ 

( , 0n,'0m) = (/ + / ) — [cos(n — m)0cos 2?7 — (— l) n sin(n — m)0sin 2r/] 

J-b Ja 2/ 

(4.11) 
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Consider the first term in the bracket [ ]. Since cos <f> is an even function of A, 
cos(n — m)(j) is even and cos2r/ is odd (see Eq. ( |2.11| )). Thus the integrand 
is odd and gives zero upon integrating over (—6, —a) U (a, b). Similarly 
the second term in the [ ] is also odd, because sin <j> is odd which implies 



sin(ra — m)4> is odd and sin 27/ is even from Eq. (2.11). Thus (ip n ,ip m ) = 
for n — m odd. 

This completes the proof of the orthonormality of ip n (X). 



5 Kernel 

We will obtain an expression for K(/i, v) analogous to the Christoffel-Darboux 
formula from Eq. (|2.5| ), which is more convenient for the present two cut 
structure. Multiplying and dividing Eq. ( |2.5| ) by fi 2 — v 2 , the kernel for 
N = IP even is 



yV /^2n+l 0)^+1 M ~ V>2n+1 W^fon+l (») i , g ^ 
n=0 7 

Using the recurrence relationship for the even and odd orthogonal polyno- 
mials and with some algebra we get 

- 1p 2 P-l(v)lp2P+l(v)) + \ZB 2 P-lA 2 p(lp2p(lJ')' l p2P-2(v) 

- ^2P-2W%p{y))\ (5.2) 
Following a similar procedure for N = 2P + 1 odd the kernel K(p,, v) is 

1 1 



K ^'^ = (2P + 1) (/j 2 - t/ 2 ) [v^P^2P+l(^2P+l(^)^2P-l(^) 

- 1p 2 P-l(fJ-)lp2P+l(v)) + V 7 B 2 p + lA 2 p{ll)2P+2{v)'<p2p{v) 

~ 1p2p(fJ-)lp2P+2(^))] (5.3) 

The kernel in this form is more useful than the standard form for the 
symmetric two cut hermitian matrix model due to the explicit form of / and 
<f> (Eqs. and (glCj) ) in which 2</> appears rather than <f). 
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A convenient form for the kernel with two cuts as derived above is (for 
N even) 



+ 



N (fi 2 - v 2 ) 
1 1 
N (n 2 - v 2 ) 



^A N B N _x (ip N {ii)ipN~2 {v) - ^n^iA^nW))]- (5.4) 



Using the asymptotic ansatz and doing some simple trigonometry with 



/Hp) 



cos(N( + (/)+x-ri)(n), ipN-i{v) 



7|=cos(iVC-0+x- 



cos(N( + x + >7)(a0 an< i i>N-2{v) 



i 



cos(iV"C 



2(j> + x + v)^) we S e t the following 

— 1 ^A N B N+ i 



where Nh^>(jj) 



[ (cos20Gu) - cos 20(i/)) [cos iV(/i (o) (/u) + /i (o) (i/)) 

+ cosiV(/i (o) (^) -/i (o) H)] 

+ (sin 20(/x) - sin 20(i/)) sin JV(/i (o) (/i) + fcW (i/)) 

+ (sin20(/i) +sin20(^))siniV(/i (o) ( / u) - h {o) (u))] (5.5) 

: (JVC + X + </>-r/)(/i) and 

— 1 ^AnBn-i 



™ Ui 2 -v 2 )y/Mf{y) 

[ (cos20(/x) - cos 20(i/)) [cos iV(/i (e) (/u) + /i (e) (i/)) 

+ cosJV(/i (e) (/z) - fr (e) (*/))] 

+ (sin20(/i) -sin2</>(^))siniV(/i (e) (^) + / l (e) (^)) 

+ (sm20(/i) +sin20H)siniV(/i (e) (/i) -h {e) (v))] (5.6) 

where Nh^fjj) = (JVC + X + ^ (/•*)• Using the fact that lim^ M if (/i, i/) = 
p(/i) one finds that C'(/i) = — 7r P(/ u )- I n m ore detail, the first term in 
K(°\fi,v) and JT (e) (/^, 

f) is of order O(^) while the second and third terms 
give a term of order 0(1) which arises due to a Taylor expansion about 
v = fi + 6. The O(l) terms each give a factor ^h'^(fi) and ^h'^(fx) on 
taking the z/ — > \i limit. These terms then combine to give —('(fi). This 
determines the unknown function C(/-0 hi the asymptotic ansatz ifi n (n) in 
terms of the density of eigenvalues, and yields Eq. (2.9). 
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Eqs. (HD and (fTS) together with K = + K^ e \ and Eqs. (O- 



2.11| ) describe the full kernel function for the hermitian one matrix model 



with eigenvalues distributed in two disjoint cuts. The result is valid for any 
polynomial potential which can support two symmetric cuts. When a = 
the two cuts in the density of eigenvalues merge into one, there is a phase 
transition from a double band to a single band. The two-cut ansatz, Eq. 
(|2.7D for ijj n (X), changes to the one-cut ansatz of ref. Q Eq. (2.6) at a = 0. 
Thus the exact kernel derived above contains all the information of phase 
transitions and universality for both the fine and coarse grained correlators. 

6 Correlation Functions 

The full two-point connected correlation function for the two-cut hermitian 
matrix model is 

p c (p,u) = -[K{p,v)f = -{[K^{p,y)f + [K^{p lV )f 

+ 2K^(p,y)K^(p,y)) (6.1) 

with K®{pL,v) &ndK^(p, v) given by Eq. (|5.5| ) and Eq. (|5.6| ) respectively. 
We now discuss the correlator in various limits. 

(a) . First consider the fine grained correlator i.e. where the separation p — v 
of the order of the eigenvalue spacing m O(jt) ; 

(p - u) « O(i) = 26 (6.2) 

Then only the third terms in K(°\n, v) and K^ e \/j,, v) contribute to K(fi, v). 
Thus (for N even) 

where fi = Over the range 5, p(p) is just constant which can be scaled 
away. This extends the validity of Dyson's short distance universality to 
eigenvalues distributed in two disconnected cuts. The two point correlation 
function is 

p^,u)=p(p)p(u)[l-( S -^) 2 ) (6.4) 

with x = 2Nir5p(p). 

(b) . For 5 << O(jf) and x « 1 one gets 

p{p,v)^p\p)[N{pL-y)}\ (6.5) 
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It is evident that the exact short distance correlator has no singularity as 

V — > (JL. 

(c). In the smoothing regime, i.e., for \i — v » O(jr) we average over /j,, v 
in an interval A such that jj << A << 1. Replacing < cos 2 Nh > and 
< sin 2 Nh > by a half and < cos Nh >=< sin Nh >= 

-<[K{im, V )] 2 > 

1 fXU 1 

~ 2N 2 ir 2 (ii 2 - v 2 ) 2 sin 2cp(fi) sin 2<t>(v) 
(1 — cos 24>(n) cos 2<p{v)){l + cos <j){p) cos 4>{y)) 
sin 2(p(fi) sin 20(i/) sin </>(//) sin 0(f)] , (6-6) 

where 0(^u) = </>(/i) — 2rj((i), is the global smoothed correlator valid all over 
the cuts. It is interesting to note that even when // and v belong to the 
same segment the smoothed correlator has changed from the single band 
universal result (ref. @] Eq. (2.22)). 

This proves the universality of smoothed global two point correlators 
(as the correlator depends only on the end points of the cuts (—6, —a) and 
(a, b)) for symmetric two-cut matrix models, generalizing the universality 
for the one cut matrix models proven in Jl9| , [I]]. Note that the function x 
undetermined in our approach is absent in the universal limits (a), (b) and 
(c). 

7 Conclusion 

In this work, the orthogonal polynomial method for large N hermitian 
matrix models with arbitrary potentials that support eigenvalues distributed 
in two symmetric cuts has been developed. An asymptotic form for the or- 
thogonal polynomials ip n (X) valid for n close to N and A inside the cuts 
has been obtained. This is done by making an ansatz in terms of unknown 
functions f,4>,rj,x and £. The first three are then constrained using the 
recursion relations and orthonormality, which determine them to be univer- 
sal functions of A independent of the form of the potential except for the 
endpoints of the cuts, x, though not determined in this approach is also 
conjectured to be universal. The crucial difference between the one-cut case 
and the two-cut case appears in the behaviour of the recursion coefficients, 
which have a "period two" character, i.e., are described by two continuous 
functions in the large N limit instead of one. This modifies the functional 



+ 
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form of the functions / and 4> and introduces the new nontrivial function rj 
in the asymptotic ansatz that vanishes in the single cut case. 

The orthogonal polynomials are then used in a Christoffel-Darboux for- 
mula (a version of this formula especially suited for two symmetric cuts is 
used) to determine the eigenvalue kernel K(fA, v). The relation between this 
kernel and the eigenvalue density fixes the last unknown function £ in the 
ansatz. This is the only non universal function in the asymptotic formula for 
the orthogonal polynomials and hence also in the kernel and the correlation 
functions. 

The kernel determines global fine grained eigenvalue n-point correlation 
functions through well known formulas. The two point function is discussed 
in detail. It is shown to be universal in two limits, the local fine grained 
limit and the global smoothed limit (in the latter it depends on the potential 
only through the endpoints of the cuts). It is hoped that the above can be 
generalized to multi-cut as well as complex matrix models. 
Note added : 

After this work was completed I became aware of f20|| , where the or- 
thogonal polynomials are obtained for the case of quartic potential V{M) = 
^-M 2 + ^M 4 with two cuts by different methods. I thank Prof. E. Brezin 
for information about this work. The method and results of the present 
paper are valid more generally, i.e., for arbitrary Z2 symmetric polynomial 
potentials. 

I would like to thank the referee for pointing out the following references: 
i) |22| for a renormalization group approach to universality of smoothed 
correlators in multi-cut models, ii) (T^j for applications of two-cut models to 
QCD, and iii) [^l| for computation of smoothed Green function for multi-cut 
models. It is not clear that the connected density-density correlator derived 



from [21] is the same as the one presented here. This may be due to the 
different limiting procedures used (e.g. the iV — > 00 and Asym — > limits 
do not commute flllD. I would also like to thank G. Akemann for an email 



discussion on this point. 
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I would like to thank E. Brezin and S. Jain for discussions. 
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